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ABSTRACT: A whole genome association (WGA) study was performed to detect significant polymorphisms for meat quaUty traits in 
an F2 cross population (N = 478) that were generated with Korean native pig sires and Landrace dams in National Livestock Research 
Institute, Songwhan, Korea. The animals were genotyped using Illumina porcine 60k SNP beadchips, in which a set of 46,865 SNPs 
were available for the WGA analyses on ten carcass quality traits; hve weight, crude protein, crude lipids, crude ash, water holding 
capacity, drip loss, shear force, CIE L, CLE a and CIE b. Phenotypes were regressed on additive and dominance effects for each SNP 
using a simple linear regression model, after adjusting for sex, sire and slaughter stage as fixed effects. With the significant SNPs for 
each trait (p<0.001), a stepwise regression procedure was applied to determine the best set of SNPs with the additive and/or dominance 
effects. A total of 106 SNPs, or quantitative trait loci (QTL) were detected, and about 32 to 66% of the total phenotypic variation was 
explained by the significant SNPs for each trait. The QTL were identified in most porcine chromosomes (SSCs), in which majority of 
the QTL were detected in SSCs 1, 2, 12, 13, 14 and 16. Several QTL clusters were identified on SSCs 12, 16 and 17, and a cluster of 
QTL influencing crude protein, crude lipid, drip loss, shear force, CIE a and CIE b were located between 20 and 29 Mb of SSC12. A 
pleiotropic QTL for drip loss, CIE L and CIE b was also detected on SSC16. These QTL need to be validated in commercial pig 
populations for genetic improvement in meat quality via marker-assisted selection. (Key Words: Whole Genome Association, SNP, 
Meat Quality, Korean Native Pig, Landrace) 



INTRODUCTION 

Consumers' preference on pork has enabled to produce 
high meat in pig industry. To meet this ends, studies on 
genetics to identify genes or chromosomal regions (QTL) 
that are associated with meat quality have been performed, 
and lots of QTLs (6,347) for economically important traits 
have been reported in pig (http://www.animalgenome.org/ 
QTLdb/). 

Experimental crosses of two genetically divergent 
breeds, such as wild boar vs. European domestic breeds, or 
Meishan vs. western breeds were extensively used in QTL 
studies (Andersson et al., 1994; Rothschild et al., 1995, 
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2007; Janss et al., 1997; Bidanel and Rothschild, 2002; Hu 
et al., 2007). After the first report on genome-wide scan and 
detection of QTL for growth, length of small intestine, and 
fat deposition in pig (Andersson et al., 1994), many QTL 
studies using F2 cross populations were performed on 
growth and fat deposition (Knott et al., 1998; Paszek et al., 
1999; Rohrer, 2000; Wada et al., 2000; Bidanel et al., 2001; 
Malek et al., 2001a; Kim et al., 2005), carcass traits 
(Andersson-Eklund et al., 1998; Rohrer et al., 1998), and 
meat quality (de Koning et al., 1999, 2001; Malek et al., 
2001b). Recently, several QTL studies on purebreds for 
meat quality were also reported (Uemoto et al., 2008; Soma 
et al.,2011). 

Korean native pig (KNP) were reported with good meat 
qualities such as better redness, less cooking loss and shear 
force than Yorkshire or Landrace breeds (Cho, 2006), and 
with good amounts of amino acid compositions that were 
possibly related to flavor and palatabiUty (Hwang et al., 
2004). Kim et al. (2007) and Choi et al. (2011) conducted 
genome-wide scans using microsatellite markers in an F2 
cross population between KNP and Landrace to detect 
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Mendelian and parent-of-origin QTL for growth, body 
composition, and meat quality traits. 

Recently, with the developments of high throughput 
genotyping technologies such as next generation 
sequencing and high density porcine SNP chips, it is 
possible to genotype large amount of SNPs, e.g. with the 
porcine 64k Illumina Infinium assay, which enabled 
genome-wide association (GWA) study to detect QTL with 
high mapping resolution up to lor 2 Mb confidence region 
(Fan et al., 2009; Gorbach et al., 2009; Onteru et al., 2009). 

In this study, we performed whole genome scans to find 
significant sets of SNPs that were related to meat quality 
using the Dlumina porcine 64k SNP chips in the KNPx 
Landrace F2 cross population that was used in Kim et al. 
(2007) and Choi etal. (2011). 

MATERIALS AND METHODS 

Animals and phenotypes 

Data were collected from the QTL experimental 
population that were produced by crossing Korean Native 
(KNP) boars and Landrace (LN) sows at National Livestock 
Research Institute (NLRI), Songhwan, Korea. Ten Fi boars 
were randomly chosen for inter se mating up to six Fi sows 
to produce 38 full-sib F2 families resulting in 490 F2 
progenies from the first (N = 281), second (N = 130), and 
third (N = 79) parities. The raising and management 
practices of the animals were described in Choy et al. 
(2002a; 2002b). The F2 individuals were slaughtered at the 
age of 201 to 251 d with the average of 214.8 d at the NLRI 
slaughterhouse. The average final (live) weight before 
slaughter was 91.6 kg with a range of 45.2 to 145 kg. Data 
measurements on ten carcass quality traits (live weight, 
crude proteins, crude lipids, crude ash, water holding 
capacity (WHC), drip loss, shear force, CIE L, CIE a, and 
CIE b, were described in detail in Choi et al. (2011). 

Molecular data 

Using the Infinium HD Assay Ultra Protocol (Illumina), 
all of the 490 F2 animals were genotyped with the porcine 
64k SNP beadchips (Illumina Inc., USA). Raw data were 
visualized and analyzed with the Genome Studio software 
(Illumina). For subsequent data analysis, a subset of SNP 
was selected by removing the SNPs in sex chromosomes. 
Every SNP from the chip data was screened for the 
availability of GWA tests. Those SNPs on autosomal 
chromosomes were removed before WGA testing, which 
met the following three criteria; i) the number of genotype 
group with one or none (e.g. only AA genotypes and no AB 
or BB), ii) with a minor allele frequency less than 0.05, and 
iii) with the proportion of genotyped individuals less than 
90%. 



Statistical analysis 

A SAS GLM procedure (SAS version 9.1) was used to 
preadjust the animal phenotypes before WGA testing. Sex, 
sire and slaughter stage were fitted as fixed effects for all 
meat quality traits, except the live weight, for which the 
effect of slaughter stage was left out of the model. Then, the 
residuals of each phenotype were regressed on the additive 
and dominance SNP effect for each SNP imder a simple 
linear regression. In the model, SNP genotypes with AA, 
AB and BB were assigned as 1,0, and -1 for the additive 
effect, and 0, 1, 0 for dominance effect, respectively. For 
significance threshold, 0.1% point- wise p value from F 
distributions was applied for each SNP test. 

Among the significant SNPs, the best set of SNP 
markers for each trait were selected using the stepwise 
regression procedures (Neter et al., 1990), because some of 
the significant SNPs would yield redundant information due 
to linkage disequilibrium (LD) between closely linked 
SNPs, i.e. a non-random association between alleles of 
different SNPs. Inclusion and exclusion of each SNP out of 
the model was determined at p<0.001 level. 

With the additive and/or dominance effects of the 
significant SNPs that were determined by the stepwise 
regression procedures, variation explained by each (/''') SNP 
iSjfjpj) was calculated as a-f^-^i^, where i denotes 

each genotype, at is the estimated additive and/or 
dominance effects that were included in the multiple SNPs 
model, e.g. —a , d , and +a for BB, AB, and AA, 
respectively, for the SNP both with additive ( a ) and 
dominance (d ) effects, f is the frequency of i"' genotype, /U 
is the population mean that can be expressed as a (/aa'/bb) 
+ df^ (Falconer and Mackay, 1996). Proportion of 
phenotypic variance due to the / SNP was then estimated 
as Sg,^pj is],, in which S] was phenotype variance that 

was obtained from residual values of the ttait after adjusting 
fixed and/or covariate effects. The proportion of phenotype 
variance due to all of the significant SNPs was estimated as 

(iRSS„„sNp-RSSsNPs)IRSS„„sNp), where RSS is residual sum of 
squares for the model with or without fitting the significant 
SNPs of each trait. 

RESULTS AND DISCUSSION 

Summary statistics for the ten meat quality traits were 
displayed in Table 1. The coefficient of variation (CV) for 
crude lipid was the greatest (96%), followed by for drip loss 
(87%). Among the color scores, CIE b had the greatest CV 
(43%), while the lowest CV was observed for crude protein 
(5%). 

A set of 46,865 SNPs was chosen from the 62,163 SNPs 
in the Illumina porcine beadchip (Table 2). The number of 
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Table 1. Summary statistics for 490 F2 plienotypes on meat quality in a KNPxLN cross population 




1 30 1 


Trait 


Average 


Std. Dev." 


Minimum 


Maximum 


CV* 


Live weight (kg) 


91.6 


15.8 (14.0) 


45.2 


145.0 


17.3 


Crude protein (%) 


22.1 


1.1 (1.0) 


16.4 


25.3 


5.0 


Crude lipid (yc) 


2.1 


2.U (l.o) 


0.1 


23.6 


96.2 


Crude ash (%) 


1.0 


0.1 (0.1) 


0.8 


2.9 


12.5 


Water holding capacity (%) 


57.5 


5.3 (4.8) 


45.9 


80.2 


9.3 


Drip loss (%) 


2.3 


2.0 (1.8) 


0.0 


11.1 


86.7 


Shear force (kg) 


3.7 


1.2(1.2) 


1.3 


9.9 


33.4 


CIEL 


50.7 


4.5 (4.0) 


34.7 


67.1 


8.8 


CIEa 


9.7 


2.5 (2.3) 


4.5 


18.3 


26.0 


CIEb 


5.1 


2.2(1.9) 


0.7 


14.3 


42.8 



' Standard deviation. The values in parenthesis indicate residual standard deviations that were obtained for the phenotypes after preadusting for 

appropriate fixed effects. 
'' Coefficient of variation (%). 



SNPs (6,622) was the greatest for Sus scrofa chromosome 
(SSC) 1, while SSC12 had the smallest number of SNPs 
(1,098). The SSCs 4, 7, 13, and 14 had more than 3,000 
SNPs. The physical map with all of the available SNPs 
spanned about 2,424 Mb with an average distance of 
44.6±133.5 Kb between adjacent SNPs. However, the 
average distances were different between chromosomes, 
ranging between 35.6 Kb in SSC17 and 87.3 Kb in SSC6. 
Only 75% of the SNPs that were embedded in the porcine 
chip were available in this study, partly due to 
ascertainment bias, i.e. the lllumina porcine 64k SNP 
beadchip was based on the genetics of other commercial 



breeds, not of Korean native pigs. 

A total of 106 SNPs were detected in various 
chromosomes by the stepwise regression procedures, and 
the sets of SNPs for the ten carcass quality traits explained 
significant proportion of phenotypic variance, ranging 
between 31.8% for WHC and 65.6% for drip loss (Table 3). 
For live weight, twelve SNPs (QTL) were found on SSCs 1, 
2, 5, 7, 13, 14, 15, and 16, of which multiple QTL were 
located on SSCs 5, 7, 13, and 14 (Table 3). Murani et al. 
(2005) performed association analyses in Pietrain and 
DurocxPietrain populations, and found QTL for ham weight 
on SSC2 (EPOR gene at 52 Mb), which was closely located 



Table 2. The numbers of available SNPs in the Korean native pig and Landrace F2 population and average interval distances between 
adjacent SNPs in all 18 Sus scrofa autosomes (SSC) 



SSC 


Number of SNPs' 


Average interval (kb) 


Standard deviation (kb) 


Total distance'' (bp) 


1 


6,622 


44.6 


133.5 


295,554,054 


2 


2,677 


57.3 


169.2 


153,324,765 


3 


2,178 


67.1 


201.1 


146,204,919 


4 


3,693 


39.1 


90.2 


144,452,577 


5 


2,257 


45.9 


141.2 


103,518,423 


6 


1,967 


87.3 


236.3 


171,673,133 


7 


3,554 


37.5 


34.9 


133,307,931 


8 


2,224 


67.5 


168.1 


149,995,098 


9 


2,638 


58.5 


189.8 


154,365,393 


10 


1,363 


57.0 


182.5 


77,744,768 


11 


1,862 


45.2 


111.3 


84,223,424 


12 


1,098 


62.2 


170.1 


68,252,608 


13 


3,633 


59.9 


168.9 


217,637,564 


14 


4,056 


36.7 


28.6 


148,678,088 


15 


2,659 


63.7 


204.8 


169,455,966 


16 


1,574 


55.5 


168.4 


87,281,785 


17 


1,672 


35.6 


56.9 


59,523,835 


18 


1,138 


51.8 


129.2 


58,905,756 




Total: 46,865 


Average: 44.6 


Average: 133.5 


Total: 295,554,054 



" Among the 62,163 SNPs in the lllumina Porcine 64k SNP beadchip, each SNP was selected for whole genome association (WGA) tests if the number of 

genotype group is two or three, the minor allele frequency was greater than 0.05, and proportion of genotyped individuals was greater than 90%. 
' The distances between the first and the last SNPs that were located on their respective chromosomes. 
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Table 3. Identities, positions of tlie SNPs with significant effects on meat quality traits in a Korean native pig and Landrace cross 
population 



Trait/SNP marker" SNP" SSC Kb' Additive^ SE' Dominance^ SE' -logioP« % 



Live weight (kg) 


















36.70 


ALGA0002905 


[A/G] 


1 


35,240 


1.42 


0.85 


4.81 


1.17 


2.19 


3.40 


DRGA0003167 


[T/C] 


2 


47,367 






4.47 


1.15 


4.00 


2.52 


ALGA0030609 


[A/G] 


5 


8,998 


2.28 


0.83 






4.16 


1.33 


ASGA0095oo4 


[A/G] 


5 


13,541 






5.81 


1.13 


4.13 


4.14 


MIGAOO 10442 


[T/G] 


7 


62,213 






5.15 


1.12 


4.14 


3.38 


MIGAOO 10269 


[T/G] 


7 


50,229 


-3.39 


0.84 






4.20 


2.66 


INRA0040139 


[T/C] 


13 


29,753 






4.79 


1.12 


4.61 


2.90 


CASI0007467 


[T/C] 


13 


95,955 






5.03 


1.29 


4.97 


2.76 


INRA0043228 


[A/G] 


14 


25,540 


2.51 


0.97 


4.46 


1.36 


3.96 


3.79 


DRGA0014720 


[T/C] 


14 


136,395 


-2.79 


0.81 


2.59 


1.20 


3.30 


3.02 


ASGA0068956 


[A/C] 


15 


14,116 


-1.39 


0.77 


4.26 


1.23 


3.44 


2.69 


H3GA0045861 


[T/C] 


16 


2,073 






4.26 


1.17 


4.36 


2.28 


Crude protein (%) 


















40.90 


INRA0000425 


[A/G] 


1 


6,735 


0.24 


0.08 


0.36 


0.10 


3.26 


5.31 


ASGA0006575 


[T/C] 


1 


221,913 


0.19 


0.07 


0.25 


0.10 


4.66 


2.69 


A C\ A f\f\f\(\ A ^ C\ 

ASGA0009469 


[T/C] 


2 


8,808 


-0.40 


0.09 


-0.52 


0.12 


2.96 


9.92 


MARC0040826 


[T/C] 


9 


37,768 






-0.46 


0.08 


5.70 


4.88 


MARC0073944 


[A/G] 


10 


37,348 






-0.35 


0.08 


3.57 


2.92 


ASGA0054479 


[T/C] 


12 


20,575 


-0.22 


0.06 






217 


2.16 


H3GA0034945 


[T/C] 


12 


28,254 


0.29 


0.06 


0.23 


0.09 


2.06 


5.19 


ALGA0067249 


[A/C] 


12 


28,672 


-0.62 


0.14 


0.39 


0.15 


7.07 


12.70 


MARC0009546 


[T/C] 


12 


65,388 


0.46 


0.06 






2.41 


10.70 


DRGA0016110 


[A/G] 


16 


20,494 


0.24 


0.05 






4.16 


2.99 


H3GA0046813 


[A/G] 


16 


36,848 


-0.14 


0.06 


0.28 


0.08 


2.26 


2.86 


Crude lipid (%) 


















57.30 


ASGA0098229 


[A/C] 


6 


31,333 


0.60 


0.15 


-0.73 


0.19 


5.43 


6.97 


MARC0057020 


[T/C] 


6 


39,767 


-0.99 


0.25 


-0.81 


0.28 


5.50 


9.54 


H3GA0022262 


[A/G] 


7 


94,000 


0.18 


0.10 


0.72 


0.14 


2.00 


4.02 


DIAS0002084 


[T/C] 


12 


14,931 


0.30 


0.14 






3.32 


0.52 


MIGAOO 16777 


[A/C] 


12 


21,752 


0.36 


0.09 






4.06 


2.12 


MARC0004712 


[A/G] 


12 


27,642 


0.46 


0.11 






2.90 


3.18 


ALGA0067216 


[T/C] 


12 


28,517 


-1.27 


0.26 


1.03 


0.32 


4.03 


10.30 


ALGA0067249 


[A/C] 


12 


28,672 


0.55 


0.24 


-0.72 


0.26 


7.07 


5.23 


ASGA0099848 


[A/G] 


12 


61,692 


-0.44 


0.12 


-0.48 


0.14 


11.04 


4.42 


MARC0037347 


[A/C] 


13 


19,655 


-1.26 


0.23 


-1.15 


0.27 


3.51 


15.50 


Crude ash (%) 


















54.50 


ALGA0008433 


[A/G] 


1 


205,467 


-0.01 


0.01 


-0.03 


0.01 


6.46 


5.60 


DRGA0002289 


[A/G] 


1 


221,787 


-0.16 


0.02 


-0.16 


0.03 


11.32 


13.1 








lUo, / jD 






-u.u^ 


ft ni 




J.jL 


ALGA0037714 


[T/C] 


6 


43,718 


-0.12 


0.02 


-0.13 


0.02 


8.07 


8.05 


ASGA0035362 


[A/G] 


7 


100,680 


-0.05 


0.01 


-0.07 


0.02 


2.76 


4.00 


M1GA0015451 


[T/C] 


11 


62,440 


-0.04 


0.01 


-0.04 


0.01 


4.84 


6.65 


MARC0055183 


[T/G] 


14 


84,908 






-0.17 


0.03 


6.32 


6.78 


ALGA0093831 


[A/G] 


17 


20,677 


0.07 


0.02 






5.71 


7.21 


MARCO 11 5266 


[T/G] 


17 


25,587 


0.18 


0.02 


-0.12 


0.02 


10.98 


10.20 


ASGA0076045 


[T/C] 


17 


25,646 


0.05 


0.01 


-0.06 


0.02 


6.31 


8.18 


H3GA0048306 


[T/C] 


17 


25,884 


-0.13 


0.02 


-0.06 


0.02 


5.82 


6.14 


H3GA0050644 


[T/C] 


18 


14,299 


0.08 


0.02 


0.06 


0.02 


2.02 


1.97 
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Table 3. Identities, positions of tlie SNPs with significant effects on meat quality traits in a Korean native pig and Landrace cross 
population (Continued) 



Trait/SNP marker" 






JVU 


AHHitivf^ 




U\J\.W\.\. 1 CUIL/C 


ope 




/o a p 


Water holding capacity (%) 


















1 1 cn 
Ji.oU 


L/KVJ/\UUUZ / 7 1 


r A /m 


z 


1 1 7A1 
11, /Ol 




n fn 
u.o / 


11^ 

-J. 1 J 


n 8^ 

U.oJ 


9 7A 

z. /o 


1 1 An 

1 J.OU 


AT Anni '^h.q 


[i/L-J 


Z 


An Q/1/1 

Du,y44 


1 <Q 

-1 .jy 


n /lo 

U.4Z 


0 ni 

Z.Ul 


n ^A 
U.JO 


A '^0 

4. jZ 


7 10 

/ . IZ 


H'^n A 009 7 ^ 




q 


98 OJ^Q 


0 9n 

u.o / 


0 90 


1 7zL 


0 A9 


4 

J J 


4 SO 


Ayf A 1? r^c\c\c\<'i 1 n 

JVI/VKL^UUUJ / 1 / 


r A ic^ 


o 
y 


jj,1j4 


1 nn 


n n 


1 1/1 

1.J4 


n /i< 

U.4j 


0 nA 
z.uo 


/I n/1 

4.U4 


u^in Ann^9/ioi 

rl Jl_rAUU JZ4y 1 


r A irx\ 


11 


J 1 SizA 


1 n< 


n iA 






< 1Q 

j.iy 


1 A7 
1.0/ 




FT /PI 




0,4jZ 






1 ^1 


n 70 

u. /y 


0 OS 

z.yo 


1 so 
j.oy 


AT Am (YTXICk 


r A /r^i 


1:? 


1 n n^n 


1 1 'X 


n 10 


1 0^ 


n /1/1 

U.44 


0 ^/l 
Z.j4 


A on 

4.ZU 


A T A nn/^Q 1 no 




IJ 


1 / ,u /o 






1 A/1 
-1 .04 


n /in 

U.4U 


/I cn 

4.oU 


0 QS 
Z.oo 


A cr^ A nnA^ i cq 
AoLjAUUoJ loo 


r A /r^i 


1 A 
14 


40,4U4 


1 "59 
l.JZ 


n OQ 
U.zy 






1 io 

j.jy 


1 oo 
j.yz 


AT r"' AnnQTC/n 


r A /r^i 


lo 


1 O QQA 

iy,yyD 


-l.jl 


n iQ 
U.Jo 


0 A/1 
-Z.04 


n /lo 
u.4y 


A ^1 
4. J / 


11 in 

1 1.IU 


unp loss ^70 ) 
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to the QTL (47 JVIb) in this study. Soma et al. (2011) 
performed genome-wide scans in Duroc and detected QTL 
for carcass length and body length at 102 cM and 105 cM 
of SSC13 that were in close distance to the QTL (96 Mb) 



for Uve weight in this study (Table 3). 

Eleven significant QTL for crude proteins were detected, 
among which multiple QTL, i.e. two, four and two QTL 
were detected on SSCs 1, 12, and 16, respectively (Table 3, 
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Table 3. Identities, positions of the SNPs witli significant effects on meat quality traits in a Korean native pig and Landrace cross 
population (Continued) 
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SNP marker annotations and their positions. Alternative nucleotides. 
Estimates of additive and dominance effects of the SNP. ° Standard error 
^ Negative logarithm of base ten for the comparison-wise p-value of the test statistic against the null hypothesis of no SNP effect under the one SNP 

model. 

Proportion of phenotypic variance explained by each SNP or across all of the significant SNPs for the given trait. 



Figure 1). Especially the three of the four QTL on SSC12 
were clustered between 20 and 29 Mb. In our previous 
linkage mapping results using the same population as in this 
study, a QTL for crude protein was detected at 0 cM of 
SSCl (Choi et al., 2011). In the proximal region (6 Mb), we 
detected one QTL for the same trait (Table 3). 

Ten QTL for crude lipid were detected on SSCs 6, 7, 12, 
and 13, of which eight QTL were located on SSC6 and 
SSCl 2. The detected QTL explained significant proportion 
(57%) of phenotypic variation (Figure 1, Table 3). Steams 
et al. (2005) detected a QTL for fat % at 19 cM of SSC6 in 
a cross population of Berkshire and Duroc. We identified 
one QTL for crude lipid at 31 Mb of the same chromosome 
(Table 3). Sanchez et al. (2007) detected a QTL for fatty 
acids (C16:l, Palmitoleic) at 40 cM of SSC12 in a 
Duroc xLarge White cross, which was located distal to the 
crude lipid QTL at 29 Mb in this study. Munoz et al. (2007) 
detected five QTL for fatty acid composition in between 1 
and 34 cM, and three QTL at between 68 and 76 cM of 
SSC 12 in an Iberian and Landrace F2 cross, in which five 
QTL for crude lipid were detected in this study (Table 3). 
Clop et al. (2003) also reported two QTL for unsaturated 
index and linolenic acid at 34 cM of SSC6 and 31 cM of 
SSC 12, respectively, in an IberianxLandrace F2 cross 



population. In the similar regions, multiple QTL for crude 
Upid were identified in this study (Table 3). 

For crude ash, 12 QTLs were detected on eight SSCs, 
among which the four QTL on SSC17 were clustered at 21 
to 26 Mb regions (Table 3). Ovilo et al. (2002) reported 
QTL for Haematin in an IberianxLandrace F2 cross 
population at 109 cM of SSC4 and 87 cM of SSC7, 
respectively. Haematin is an iron-containing compound 
which might be related to the crude ash. We detected two 
QTL for crude ash at 109 Mb and 101 Mb of the respective 
chromosomes (Table 3). 

Ten QTL for WHC were identified on six SSCs, among 
which multiple QTL were observed on SSCs 2, 9, and 13 
(Table 3). Malek et al. (2001b) detected two QTL for the 
trait on SSC2 and SSC13 in a cross population of Berkshire 
and Yorkshire. The QTL locations were at 71 cM and 33 cM, 
respectively, that were close to the QTL in this study (Table 
3). Murani et al. (2005) detected a drip loss QTL in EPOR 
gene region (on SSC2, 52 Mb) in Pietrain, and 
Duroc xPietrain cross population. Usually drip loss and 
WHC are considered to be highly correlated each other, and 
we detected one QTL for WHC at 61 Mb of the same 
chromosome. 
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Figure 1. Test statistic profiles for crude protein, crude lipid, crude asli and Cie a under tlie additive and dominance model for eacli SNP. 
Circles indicate highly significant QTL. Y axis indicates -logio(p-value under the null hypothesis of no SNP effect). 

Nine QTL for drip loss were found on SSCs 4, 5, 9, 10, population. In the same chromosomes, two QTL for the trait 

12, 14, 15 and 16, respectively (Table 3). Jennen et al. were detected (Table 3). Li et al. (2011) reported that a SNP 

(2007) identified a candidate gene for drip loss, ATF4, at 50 in DLX3 gene at 24 Mb of SSC12 was associated with drip 

Mb of SSC5, which was close (45 Mb) to the QTL in this loss in the same population as in this study, while a QTL for 

study. Thomsen et al. (2004) found QTL for drip loss on drip loss was detected at 28 Mb of SSC12 (Table 3). In a 

SSC5 and SSC9 in a Berkshire and Yorkshire cross Duroc purebred population, Li et al. (2010) found a QTL 
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Figure 1. Test statistic profiles for crude protein, crude lipid, crude ash and Cie a under the additive and dominance model for each SNP. 
Circles indicate highly significant QTL (Continued). Y axis indicates -logio(p-value under the null hypothesis of no SNP effect). 



for the trait on SSC15 (63 cM), and Murani et al. (2005) 
also reported that the TTN gene (at 79 Mb, Ensembl 
database) on SSC15 was associated with drip loss in a 
DurocxPietrain cross population. In this study, a QTL for 
drip loss was detected at 71 Mb (Table 3). 



Eleven QTLs for shear force were detected on eight 
SSCs, of which multiple QTL were located on SSCs 1, 3 
and 13 (Table 3). Choi et al. (2011) reported a QTL for 
shear force at 19 cM of SSC15 that was detected under a 
line-cross model based linkage analyses using the same 
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population as in this study. Thomsen et al. (2004) also 
detected one QTL for tenderness at 13 cM of the same SSC 
in a BerkshirexYorkshire cross population. In this study, 
one QTL for shear force was detected at 9 Mb of the same 
chromosome (Table 3). Murani et al. (2005) reported the 
PDGFRA gene on SSC8 (35 Mb) in a Duroc cross 
population, and the NMEl gene on SSC12 (25 Mb) in 
Duroc xPietrain crosses, which were associated with shear 
force. Results from Li et al. (2011) and Choi et al. (2011), in 
which the same KNPxLN cross population was used as in 
this study, included two QTL for shear force on SSC12 and 
SSC15, respectively, which were closely located to the two 
QTL (23 Mb, 9 Mb) on the respective SSCs in this study 
(Table 3). 

Most QTL for meat color traits (CIE L, a, b) were found 
in many chromosomes, especially on SSCs 1, 2, 4, 5, 11, 12, 
14, and 16 (Table 3). In some SSC regions, many QTL for 
meat color were reported; on SSC4 and SSC7 for Minolta a 
(Ovilo et al., 2002), on SSC2 and SSC14 for Hunter L, and 
on SSC5 for Hunter b (Roher et al., 2005), on SSC2 for A 
value (Streans et al., 2005), on SSC4 for color 24hSM (Ma 
et al., 2009), on SSCl for L value, on SSC14 for A value, 
and on SSC4 and SSCS for B value (Li et al., 2010), and on 
SSC14 for CIE L (Li et al., 2011). Our study also detected 
QTL for CIE a, b, and L in these QTL regions (Table 3). 
Duthie et al. (2008) reported an epistatic effect of OPTO- 
STAR for CIE L and b that was located on SSC4, which 
also resided in the similar region to the QTL for CIE L and 
b in this study (Table 3). The QTL for CIE a that was 
located at 27 Mb of SSCl 2 had the greatest statistical 
evidence of this QTL study (Table 3, Figure 1). 

On SSCI2, many QTL were clustered between 20 and 
29 Mb for crude protein, crude lipid, drip loss, shear force, 
CIE a and CIE b (Table 3), suggesting that more than one 
gene reside in the regions. Also, on SSC16 and SSC17, 
SNPs with pleiotrophic effects were detected, e.g. 
ALGA0090076 on SSC16 for drip loss, CIE L, and CIE b 
(Table 3). 

Previously, we performed genome-scans to detect QTL 
for meat quality using the same KNPxLN cross population 
as in this study (Choi et al., 2011; Li et al., 2011). However, 
very limited QTL were confirmed in this study, partly 
because different QTL mapping approaches were apphed 
with the use of different set of markers. In this study, WGA 
was applied using high density SNP chips, while Choi et al. 
(2011) performed QTL analysis with limited numbers of 
microsatellite markers under the interval mapping models 
that were based on breed-cross designs with Mendelian and 
non-Mendelian inheritance mode of the QTL. Also, Li et al. 
(2011) performed association analyses using limited 
number of SNPs of candidate genes. 



CONCLUSION 

Selection for high growth and lean meat has been 
extensively implemented in pig industry. However, recent 
genetic improvement has been focused on better meat 
qualities. In this study, lots of QTL (108 SNPs) were 
detected for meat quality, providing a set of genetic markers 
to implement marker-assisted selection (MAS). Several 
QTL clusters were identified on SSCs 12, 16, and 17, as 
well as QTL with pleiotrophic effects, e.g. the QTL for drip 
loss, CIE L and CIE b on SSC 16. Characterization of these 
QTLs is needed to find causal mutations around the QTL, 
which would enable MAS to be implemented with greater 
efficiency and accuracy. 
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